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We study the quantum propagator in the semiclassical limit with sharp confining potentials. 
Including the energy-dependent scattering phase due to sharp confining potential, the modified Van 
Vleck's formula is derived. We also discuss the close relations among quantum statistics, discrete 
gauge symmetry, and hard-wall constraints. Most of all, we formulate a new quantization rule 
that applies to both smooth and sharp boundary potentials. It provides an easy way to compute 
quantized energies in the semiclassical limit and is extremely useful for many physical systems. 
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I. INTRODUCTION 



The most straightforward method to obtain the bound- 
state spectrum for a quantum system is to solve the 
Schrodinger equation. However, if the potential profile 
is smooth compared with the wavelength of the particle, 
the energy spectrum can be obtained by the semiclas- 
sical WKB approximation. The semiclassical approach 
reduces the task of solving the differential equation into 
a simple integral. While the great simplification is attrac- 
tive, it does not work all the times. In quantum dots, the 
confining potential is usually sharp and leads to strong 
quantum interferences which invalidate the semiclassical 
approach. One notices that the semiclassical approxima- 
tion only breaks down near the sharp confining potential. 
It motivates us to generalize the conventional semiclas- 
sical approach by including the quantum interferences 
exactly near the turning points where the semiclassical 
approximation is not appropriate. Rather nicely, we were 
able to capture the complicated quantum interference ef- 
fects by a simple energy-dependent scattering phase cor- 
rection. 

To elucidate this point, it is convenient to adapt the 
path integral formulism. Path integral provides an alter- 
native approach to formulate quantum mechanics P, 0- 
The quantum propagator G{x,x' ';T), that is the key 
quantity in quantum mechanics, is shown to equal the 
summation over all possible paths with the same end 
points. In the semiclassical limit, the dominant contri- 
bution comes from classical trajectories and fluctuations 
around them|3|, U, |5j- Within the stationary phase ap- 
proximation including fluctuations up to the quadratic 
order, the quantum propagator can be approximated by 
the Van Vleck's formula|6|. In general, there would be 
many classical trajectories that satisfies the same bound- 
ary conditions, and the Berry phase interferences be- 



tween them are important |3, la 13- By Morse's theorem, 
the second variation, considered as quadratic fluctuations 
around a given trajectory from x' to x in time T, has as 
many negative eigenvalues as there are conjugate (turn- 
ing) points along the trajectory. These conjugate points 
give rise to a Berry phase vir/2 for the trajectory, where 
v is the total number of conjugate points along the tra- 
jectory, or sometimes referred to as the Maslov or Morse 
index[l(| ■ 

Not only elucidating the crossover between classical 
and quantum mechanics, the semiclassical limit also pro- 
vides a convenient way to calculate the bound state en- 
ergy. Instead of solving the Schrodinger equation di- 
rectly, the bound state spectrum can also be computed by 
the WKB approximation ^l|- I n order to account for the 
interference effects among classical trajectories correctly, 
we rederive Van Vleck's formula with an extra scattering 
phase correction due to sharp confining potentials. 

Following the standard stationary phase approxima- 
tion and making a Legendre transformation of the time 
variable in the quantum propagator to the energy vari- 
able, we are able to generalize the Einstein-Brillouin- 
Keller (EBK) quantization rule [12, [H Q with an ad- 
ditional phase correction term 



y/2m\E - V(x)]dx = 2mr + ^ <t>s(E), (1) 



where 4> S {E) is the energy-dependent scattering phase 
due to collisions with the confining potential. The usual 
WKB approximation is the special case where the scat- 
tering phase at each turning point is assume to take on 
the energy-independent value 4> S {E) — ir/2. On the other 
hand, if the confining potential becomes infinitely sharp 
(hard- wall limit), the scattering phase rises to tt. The 
modified EBK quantization rule in Eq. (JJJ relaxes the 
requirement of the potential smoothness in the WKB 
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approximation. This is of great advantage because many 
physical systems including quantum dots, quantum wells, 
Hall bars, electronic wave guides, etc., have both hard- 
wall-like potentials (from sample edges) as well as smooth 
potentials (by applying external fields) at the same time. 
The paper is organized in the following way. In Sec. 

II, we introduce the Van Vleck's formula and apply it to 
simple systems. We explicitly show that the Van Vleck's 
approximation is incorrect in the presence of hard walls 
and the scattering phase correction is crucial. In Sec 

III, we compute the energy dependence of the scattering 
phase at each turning point and derive the modified Van 
Vleck's formula. In Sec. IV, we derive the key result of 
this paper - the modified EBK quantization rule. We 
apply it to physical systems with both smooth and hard 
confinement potentials and show that the modified term 
is necessary to obtain the correct energy levels. Finally, 
in Sec V, we relate the connection between the quan- 
tum statistics, discrete gauge symmetry to the scattering 
phase approach. Then a brief conclusion follows. 



II. QUANTUM PROPAGATOR AND 
CLASSICAL TRAJECTORIES 

In the path integral formalism^]], the quantum prop- 
agator equals the sum over all possible paths with the 
same end points, 

G(x,x';T) = (x\e~ lHT \x'} 

X>[a;]exp (i [ L(x,x,t)dt] , (2) 



where the measure T>[x] denotes all possible paths with 
end points x(0) = x' and x(T) = x. In the semiclas- 
sical limit, the phase inside the path integral oscillates 
rapidly except in the neighborhood of the classical tra- 
jectories. Within the stationary phase approximation in- 
cluding fluctuations up to quadratic order, the propaga- 
tor is approximated by the Van Vleck's formula, 



G(x,x';T) 




C p exp[iv4 p 



(3) 



where A p (x, x'; T) is the action of the classical trajectory 
starting from a;(0) = x' and ending at x(T) — x, and the 
subscript p denotes all classical paths with the desired 
end points. The strength of the quadratic fluctuations 
around the classical trajectory is 



Cp — 



d 2 A 



dxdx' 



(4) 



Finally, the total number of conjugate (or turning) points 
along the classical trajectory is denoted by v. Notice 
that, for each conjugate point, there is a 7r/2 Berry phase 
associated with it. Van Vleck's formula provides a com- 
pletely classical approximation of the quantum propaga- 
tor, in the sense that all relevant elements can be com- 
puted from the classical trajectories. 



r 



FIG. 1: Classical trajectories in the presence of two hard 
walls. On the left is a trajectory with even reflection points 
r = 2, while the right with odd reflection points r — 3. 



A straightforward example of the Van Vleck's formula 
is a free particle moving on the a finite ring with length 
L. There are infinite classical paths which satisfy the 
conditions x(0) = x' and x(T) = x. The total (route) 
distance of each classical trajectory ' + nL, 

where n is an integer. The action for each trajectory is 



A n (x, x'; T) — yrp{ x ' 



nL) 



(5) 



Taking the derivative of the action, the strength of fluctu- 
ations around each trajectory C„ = m/T is independent 
of the end points and the choice of trajectories. Since 
the particle moves at constant velocity, it is obvious that 
there is no conjugate point along any classical trajec- 
tory and thus v n — 0. Besides, because the fluctuations 
of the classical trajectory of a free particle are exactly 
quadratic, we expect the Van Vleck's formula to be ex- 
act for this system, 



G(x, x; T) 



exp 



m 
i^r^yx - 



2T V 



+ nLf (6) 



This infinite sum can be re-written in terms of its Fourier 
function with the use of Poisson summation formula in 
Appendix A.^f| Notice that 



f(y) = e la( -y + ^ 2 <-> Ftp) = J ™ e ~^/i a+l kp (?) 

V a 

Choosing a — L, the summation over coordinate y = na 
can be turned into the summation over momentum k n — 
2nir/L. The propagator is then 

G(x, x; T) = — exp[ik n (x — x 1 ) — iE n T], (8) 



where k n = 2nir/L is the quantized momentum and 
E n = k^/2m is the quantized energy. It is obvious that 
the propagator G(x, x'; T) calculated by the Van Vleck's 
formula is exact in this case. 

Let us now apply the Van Vleck's formula to another 
physical system - a free particle bouncing back and forth 
between two hard walls. We calculate the propagator 
explicitly and show that the Van Vleck's formula leads 
to incorrect results. 
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The trajectories in this problem can be classified by 
the number of collisions with the hard walls, as seen in 
Fig. 1. For those trajectories that collide with the hard 
walls even times, the route distance is d e n — x — x 1 + 2nL, 
while the distance is d° = x + x' + 2nL for trajectories 
that collide with the walls odd times. The action for each 
trajectory can be computed straightforwardly 



The Hamiltonian is 



A^(x, x'; T) 



2T 



(x-x' + 2nL) 2 , 



A° n (x,x';T) = —{x- 



2nL) 2 



(9) 
(10) 



Here A e /°(x, x'\ T) denotes the action for trajectories 
with even/odd reflection points. The fluctuations along 
all trajectories contribute the same C n = m/T as in the 
previous example. For an one-dimensional motion, a con- 
jugate point is identified as the position where the veloc- 
ity vanishes. However, for a free particle bouncing back 
and forth between two hard walls, the velocity is constant 
up to a minus sign and does not vanish at any point along 
the classical trajectory. Thus, the number of conjugate 
points is zero, v — 0. 

The propagator without any phase correction is 



Gvv(x,x';T) 



2niT 



x-x' + 2nVf 



+ exp[i — {x + x' + 2nL) 2 ] 



(11) 



Both infinite sums can be turned into summations over 
discrete momentum again by mean of the Poisson sum- 
mation formula. The prefactors cancel as in the previous 
example and we are left with the simple result, 



1 °° 

G vv (x,x';T) = — E exp[-iE n t}: 



71=0 



x I cos[k n (x — x')] + cos[k n (x + (12) 

where k n = nn/L is the quantized momentum and E n = 
k 2 /2m is the quantized energy. Combining two cosines 
would leads to cos(k n x) cos(k n x'), while the correct form 
should be sin(k n x) sm(k n x'). In fact, one can recover the 
exact answer (with all prefactors right!) if we change 
the sign of the second term in Eq. I|12|). That is, only 
if we assign an extra phase n to trajectories with odd 
reflection points, will the modified Van Vleck's formula 
become correct! 

In the following section, we study the path integral 
formalism in the presence of a single hard-wall boundary 
and show that an extra phase correction arises naturally 
due to collisions with the confining potential. 



III. SCATTERING PHASE DUE TO HARD 
WALL 

Consider a particle moving under the influence of a 
regular potential V(x) and a hard-wall potential V c (x). 



H = ^— + V(x) + V c (x), 
2m 

where V c (x) is the hard- wall potential at x — 0, 



V c (x) 



0, x > 0; 
co, x < 0. 



(13) 



(14) 



The regular potential is treated in the ordinary way while 
the hard-wall one is viewed as the depletion of Hilbcrt 
space. The complete set of the Hilbert space is now re- 
duced, 



dr\r)(r\ = 1, 



dp 



E / g^bXe^l = 1. 



(15) 
(16) 



It would become clear later that the phase <f> is associated 
with the scattering phase in the path integral. Slicing 
the time interval T into infinitesimal pieces and inserting 
complete sets of the coordinate space, the propagator is 

G(r,r';T) = (r\ e - lHT \r') 

r oo N-l 

= / dr n T[(r n+1 \e- UH \r n ), (17) 

J n=0 

where — r and ro = r' are all positive. Each ma- 
trix element in the product is computed by inserting the 
complete set in momentum space into Ea. lfTf 



{r n +i\ e 



cH 



f dp n r ■ u 1 

/ ~2~w 6XP t~ l€ n * 



E 



> e 



ip n (r n + 1 -x n )-i4 



(18) 



where the phase = for x n = r n , and <j) — n when 
x n = —r n . Since x n = ±r„, the two terms can be com- 
bined and lead to the unconstraint integral over x n . After 
changing the constrained variable r n to x n , it is conve- 
nient to write the Berry phase <p n in the following way 



4> n = 7r[6(a: n+ i) - 9(x n )]. 



(19) 



Notice that the Berry phase is zero if the path does not 
pass through x = in the infinitesimal time interval dt n 
and 7r if the path passes through. The integral over mo- 
mentum can be carried out easily and the propagator is 



G(r,r';T) 



E ^ 

x' — ±r' 



T>[x]exp[iA(r,x';T)]. (20) 



The total phase </> s = 7r[6(r) — Q(x')) is a boundary term 
and can be pulled out of the path integralQ. The paths 
are divided into two topologically distinct classes. For 
all possible paths starting from r to r', the scattering 
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phase is zero, while for those starting from r to — r', the 
scattering phase is ir that causes a minus sign. 

The classical trajectories among the paths can be then 
classified in the same way. Furthermore, trajectories with 
end points r and r' can be identified as trajectories (in 
the physical half plane) with even reflection points and 
those with end points r and — r' are trajectories with odd 
reflection points. Therefore, in the semiclassical limit, 
the Van Vleck's formula is modified with an extra phase 
term, 

G(r,r';T)~^L=VV^exp[a p -%], (21) 

V27TI 

p 

The proof for more than one turning point is straight- 
forward and the scattering phase just add up. It would 
become clear in the following section that the scattering 
phase correction is crucially important in determine the 
energy spectrum. 

IV. MODIFIED EBK QUANTIZATION RULE 

The most powerful use of Van Vleck's formula is that 
it leads to the EBK quantization rule in the semiclas- 
sical limit. One notices that, if we set x — x' in the 
propagator and integrate over all possible x, it results in 
the quantum partition function Z(T) = exp[— iE n T]. 
The energy levels can then be identified as the singular- 
ities of Z(u>) which is the Fourier transformation of the 
partition function. Within stationary phase approxima- 
tion, it can be shown that the total phase <f pdq — ivir/2 
(in the absence of sharp boundaries) is quantized|2j and 
leads to the EBK quantization rule, 

<j) pdq = 2mr + u^, (22) 

where v is the number of turning points along the peri- 
odic orbit. The usual WKB approximation is the special 
case with two conjugate points v = 2. The presence 
of the sharp boundaries changes the scattering phase at 
each turning point from n/2 to n and leads to the mod- 
ified EBK quantization rule. It is interesting to see that 
the scattering due to sharp confining potential modified 
the spectrum only through the scattering phase <j> a . 

Now we are ready to consider the confining potential 
in more general form 

V c (x) = G(-x)[V + Vi\x\], (23) 

where Vq = k\/2m is the potential height and V\ — 
k\/2m is the slope of the confining potential. The scat- 
tering due to V c (x) can be solved exactly and the eigen- 
states are 

\i/>(k)) =\k)+e- i *>W\-k). (24) 

The scattering phase is apparently energy-dependent as 
shown in Figure 2. For the hard-wall potential (k/ko = 




FIG. 2: Scattering phase for different potential height Vb = 
kg /2m and slop Vi = k\/2m. 

0, k/ki = 0), the scattering phase is n, while for the 
smooth potential (k/ko ^ 1, fe/ifei 3> 1) the phase be- 
comes 7r/2 as in the WKB approximation. Following 
similar calculation in previous section, we arrive at the 
modified EBK quantization integral in Eq. ^ 

We apply the modified EBK quantization rule to a 
finite potential well of length L and with height Vb = 
k^ /2m. After some algebra, the scattering phase is shown 
to be 4>(k) — 2 cos" 1 [(k/ko) 2 — !]• The quantized energy 
E n = k 2 l /2m satisfies 

2k n L = 2mr + (/>(k n ). (25) 

Quite surprisingly, the spectrum obtained by the semi- 
classical approach is identical to the exact solution. This 
shows that the quantum interference effects arose from 
the sharp confining potential can be captured by the scat- 
tering phase rather well. 

The modified EBK quantization rule can also be ap- 
plied to physical systems in higher dimensions. Let us 
consider a spherical or hemispherical qunatum dot. We 
can either apply the modified EBK formula directly to 
the true three-dimensional traiectories|l8j or apply the 
formula after reducing the system to one dimension. 
Here we adapt the second approach. After separation of 
variables, the radial effective Hamiltonian of the three- 
dimensional spherical (hemispherical) quantum dot be- 
comes one-dimensional with the effective potential 

(_ oo, r > a 

where I is the quantized angular momentum. For the 
spherical quantum dot, / takes on all integer values, while 
for the hemispherical dot, only odd integers are allowed 
due to the flat boundary. 

The classical trajectory of the electron is confined be- 
tween the hard- wall boundary at the surface and the cen- 
trifugal potential near the origin. Thus, there are one re- 
flection point 4> s = 7T and one conjugate point <p s = tt/2. 
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Applying the modified EBK quantization rule, the ap- 
proximate energy satisfies the algebraic equation, 



1 



(27) 



where r E = yUj + 1)/ (2mE) is the conjugate point 
and a is the radius of the dot. Instead of solving the 
Schrodinger equation directly, the energy levels can be 
determined easily by the algebraic equation in Eq. 1271 In 
the semiclassical limit, the conjugate point is close to the 
origin, i.e., a/r E 2> 1. The approximate expression can 
be further simplified, 



E n ,, 



2ma 2 



(28) 



where V = ^1(1 + 1). 

Notice that this problem can be solved exactly by the 
spherical Bessel functions. The hard-wall boundary re- 
quires the wave function vanishes at the surface of the 
sphere, ji (V 2mEa) = 0, that leads to quantized energy 
levels. In the same limit a/r E S> 1, the spherical Bessel 
function is approximated by the asymptotic expansion 
that leads to 



E< 



2ma 2 



I 

n+ 2 



n 2 



(29) 



The above exact result does not seem to agree with 
Eq. (|28|l at first glance. However, if the angular mo- 
mentum is also semiclassical (Z 3> 1), the last term in 
Eq. (2BJ} is V j2 ~ I j2 + 1/4 up to 0(1/1) corrections. It 
is then clear that both give the same result. We empha- 
size again that the agreement is only possible when the 
appropriate scattering phase is included. 

Another way to obtain the modified EBK quantization 
rule for the 1/r potential is the conventional Langer's 
correction approach. Instead of including the appropri- 
ate scattering phase, one can obtain the same energy 
spectrum by modified the potential appropriately. While 
both approaches give the same spectrum, it is known 
that the wave functions calculation in scattering phase 
approximation is more accurate. [lH ITtI ] 



V. MIRROR PROJECTION 

In the previous section, we treat the hard-wall bound- 
ary as depletion of the Hilbert space. An alternative way 
is to view it as a discrete Zi gauge symmetry of the wave 
function 



ip(x) 



(30) 



The minus sign is chosen here to make the wave func- 
tion vanishes at x = so that the boundary condition 
-0(0) = is always satisfied. Since the propagator can 
be written down as the summation of eigenfunctions, 





(a) 



(b) 



FIG. 3: Classical trajectories of two particles whose quantum 
statistics is replaced by the equivalent hard wall at x\ — X2 
in the configuration space. In part (a), a direct trajectory 
is shown and, in part (b), the shown reflected trajectory is 
equivalent to exchanging two particles which results in an 
extra Berry phase. 



G(x,x';T) = J2 n iPnix^nix 1 ) e: vp[-iE n T}, where ip n (x) 
is the eigenfunction with eigenenergy E n . The discrete 
gauge symmetry of the wave function implies that the 
quantum propagator has the symmetry 



G{x,x';T) = -G(x,-x';T). 



(31) 



Now choose both x — r and x' = r' to be positive, 
the propagator can also be viewed as the wave func- 
tion G(r 7 r'-T) = ijj r i(r,t) that satisfies the Schrodinger 
equation with a delta function source at (x,t) — (r',0). 
The propagator Go(r, r'; T) without the hard- wall bound- 
ary satisfies exactly the same differential equation except 
that the boundary condition at x = is not met. Notice 
that the mirrored propagator Go(r, r';T) = Gq(t, —r'; T) 
satisfies the Schrodinger equation without the source 
term since the delta function S(r + r') = for positive 
coordinates. Therefore, the propagator that satisfies the 
correct boundary condition is constructed as 



G(x,x';T) = G (x,x';T)-Go(x,x';T). 



(32) 



The above result is equivalent to Eq. (|20|l . It is obvious 
that the discrete gauge symmetry in Eq. (|31|l is satisfied. 
This method is just the familiar mirror charge trick in 
the classical electromagnetism. 

Since we can solve the hard-wall boundary by dis- 
crete gauge symmetry, we might as well go the other 
way around. It is possible to replace the quantum statis- 
tics between particles by the hard- wall boundaries. Let 
us consider the simplest case - two interacting particles 
with either bosonic or fcrmionic statistics. The discrete 
gauge redundancy is 



ip(x) 



i<t> 



ip(-x), 



(33) 



where x = x\ — xi is the relative displacement between 
two particles. The phase correction is <fr = for bosons 
and 7r for fermions. The discrete gauge symmetry is re- 
moved by imposing a hard wall x\ = x-i in the configura- 
tion space, and a Berry phase <f> accumulates upon each 
reflection due to the hard wall. 
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Classical trajectories are classified into two categories 
- the direct path and the reflected one as shown in Figure 
3. If we extend the reflected trajectory into the unphysi- 
cal regime inside the hard wall, as shown in Fig. 3(b), the 
reflected trajectory is equivalent to an exchange between 
two particles. This approach would be useful when study- 
ing few interacting quantum particles, e.g., two strongly 
interacting bosons or fermions bouncing back and forth 
between two hard walls. In the semiclassical limit, we 
can safely ignore the quantum statistics by solving all 
classical trajectories inside a specific triangle in the two- 
dimensional configuration space. 

VI. CONCLUSIONS 

In this paper, we study the scattering phase of classical 
trajectories due to sharp confining potentials. Inclusion 
of the energy-dependent scattering phase, the modified 
EBK quantization rule is derived. We also relate the 
hard wall boundary approach to the quantum statistics 
and the discrete gauge symmetry. Unlike the WKB ap- 
proximation that is only applicable to smooth potential 
profiles, the new quantization rule provides us with an 
easy way to estimate the energy levels in the presence of 
both smooth and sharp confinement potentials. 

We thank Darwin Chang for fruitful discussions, es- 
pecially on the mirror projection and the discrete gauge 
symmetry. This work was supported by the National 
Science Council of Taiwan, R.O.C.. 



sider a physical system on a finite ring with length L and 
lattice constant a. The total number of sites is N = L/a. 
The discrete version of the usual delta function is 



£^ = © E ^, (ad 

x=na * ' G=2nir/a 



where G is the reciprocal lattice vector. Consider the 
following summation, 



£/M= f ^F(k) e ikx , (A2) 

n ' x—na 



where x n = na and F(k) is the Fourier transformation 
of f(x). With the help of the identity in Ea. HAlf) . the 
summation over coordinates is turned into another sum- 
mation over reciprocal momenta. Taking the thermody- 
namical limit L — ► oo, the discrete delta functions are 
related to the continuous ones by LSk,c — 2n5(k — G). 
Finally, we arrive at the useful Poisson summation for- 
mula, 



APPENDIX A: POISSON SUMMATION 
FORMULA 



£/(na) = i5>(^). 

* J n — ' n 



(A3) 



Poisson summation formula provides a convenient way 
to related two infinite summations together. Let us con- 
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